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■ i-l , 1 INTRODUCTION 



ABSTRACT 

Equations of fully general relativistic radiation hydrodynamics around a rotating black 
hole are derived by using the Kerr-Schild coordinate where there is no coordinate 
singularity at the event horizon. Since the radiation interacts with matter moving with 
relativistic velocities near the event horizon, the interplay between the radiation and 
the matter should be described fully relativistically. In the formalism used in this study, 
while the interactions between matter and radiation are introduced in the comoving 
frame, the equations and the equations and the derivatives for the description of the 
global evolution of both matter and the radiation are given in the Kerr-Schild frame 
(KSF) which is a frame fixed to the coordinate describing the central black hole. 
As a frame fixed to the coordinate, we use the locally non-rotating reference frame 
(LNRF) representing a radially falling frame when the Kerr-Schild coordinate is used. 
Around the rotating black hole, both the matter and the radiation are affected by the 
frame-dragging effects. 

Key words: accretion: accretion disks — black hole physics — hydrodynamics — ra- 
diative transfer — relativity. 



» I ■ It is widely believed that the accretion flow onto black holes plays one of the essential roles in the active phenomena in the 
universe. Relativistic effects in accretion disks around black holes derives some major activities of astrophysical black holes, 
such as active galactic nuclei, Galactic black hole candidates and possibly gamma-ray bursts. For the mass accretion rate near 
or over the Eddington mass accretion rate, the interactions between radiation and matter in the accretion disk are important. 
While for the supercritical accretion flows considered in e.g. Syfert galaxies or the black hole binaries the photons interact 
with matter, for the hypercritical accretion flows considered in the central engine of gamma-ray bursts the neutrinos interact 
with matter around the central black hole. In such situations, the dynamics and the energy balance in matter and radiation 
are affected by each other. 

So far, the general relativistic radiative transfer is investigated by many authors (e.g. Lindquist 1966; Anderson & Spiegel 
1972; Schmid-Burgk 1978; Thorne 1981; Schinder, Bludm an & Piran 1988; Turolla & Nobili 1988; Anile & Romano 1992; 
Cardall & Mezzacappa 2003; Park 2006). Linda uist 1 1966h gives a general treatment of the radiation transfer equat ion and 



the radiation hydrodynamic equations derived by using a comoving Lagrangian frame of reference. Thorne! (1981) derive 



general relativistic moment equations up to an arbitrary order by introducing projected symmetric trace-free (PSTF) tensors. 
Since these formalisms are based on the comoving frames, the physical quantities in terms of matter and radiation and 
the directional derivatives are also described in the comoving frame. In some astrophysical objects, the radiation interacts 
with matters moving at relativistic velocities. In such cases, the interaction between the radiation and the matter should be 
described fully relativistically. This is easily done if both radiation field and mat t ers ar e evaluated in the comoving frame which 
is a frame where the element of matter is at rest. Based on this idea, Mihalasl 1 1980h introduce the radiation hydrodyr 



/namic 



equations in the Eulerian framework. In this approach, the physical quantities in terms of matter and radiation and the 
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derivatives were introduced in the frame fixed to the coordinate of the central object, e.g. a black hole, while the interactions 
between matter and radiation were calculated in the comoving frame. That is, the local processes like the interaction between 
the matter and the radiations are evaluated in the comoving frame, while the derivatives which are used for the calculations 
of the global dynamics of both the matter and the radiation are derived in the frame fixed to the coordinate describing the 
central object. 

In principle, this formalism can be extended to arb i trary space- time. IParkl (jl993r i derived the radiation hydrodynamic 
equations for spherically symmetric systems, and Park (2006) give the explicit expressions for the basic equations of the 
general relativistic radiation hydrodynamics in Schwarzs child space-time. In addition, the equations of the general relativistic 
radiation hydrodynamics in Kerr space-time are given bv lTakahashl d2007bh based on the Boyer-Lindquist coordinate which is 
one of the traditional coordinate of the metric of a rotating black hole and is most frequently used in the past studies. However, 
it is known that the Boyer-Lindquist coordinate have the coordinate singularity at the event horizon. Due to this coordinate 
singularity, the flow structure calculated by using this coordinate exhibit some unrealistic behavior near the event horizon. 
For example, the accretion flow plunges into the horizon with the speed of light, and the angular velocity of the accretion 
flow become equal to the angular velocity of the frame-dragging at the event horizon (e.g. Takahashi 2007a). Both these 
features are caused by the coordinate singularity, and the actual flows near the horizon do not have such properties. So, the 
alternative coordinate for the description of the metric around a rotating black hole are sometimes used mainly in the dynamical 
hydrodynamic or magnetohydrodynamic simulations JPapadopoulos fc Font! 19981: Font. Ibanez fe Papadopoulos 19981: Cook 



2000l ; iFonthood ; iKomissarovlboOll ; iGammie. McKinnev fc T6thl2003l ; lKomissarovH2004l ; lGammie. Shapiro fc McKinnevi2004h 
One of such coordinates is the Kerr-Schild coordinate where there is no coordinate singularity at the event horizon. By using 
this coordinate, the structure of the accretion flow can be calculated accurately near or just on the event horizon. So, since 
the inner boundary condition can be set inside the horizon, the equations of the radiation hydrodynamics around a rotating 
black hole based on the Kerr-Schild coordinate become useful when the numerical calculations are performed. One of the 
main purposes in this paper is to give these equations. 

After giving the general forms of the basic equations for the general relativistic radiation hydrodynamics in §2, in §3 by 
using the orthonormal tetrads fixed to the coordinate (§3.1), the radiation moments (§3.2), the radiation four-force (§3.3) are 
derived. In §4, we also give the basic equations for the radiation hydrodynamics including the continuity equation (§4.1), the 
hydrodynamic equations (§4.2) and the radiation moment equations (§4.3). Concluding remarks are given in the last section. 
In this paper, we assume c = 1 in most equations except in a few cases where c is explicitly used for clarity. Latin and Greek 
indices denote spatial components and spacetime components, respectively. V M denote the covariant derivative with respect 
to g^u. 



2 COVARIANT EQUATIONS FOR GENERAL RELATIVISTIC RADIATION HYDRODYNAMICS 

Before the explicit expressions for the equations of the radiation hydrodynamics in Kerr space-time are given, here we briefly 
summarize the the basic equations of the general relativistic radiation hydrodynamics in covariant form (e.g. Mihalas & 
Miharas 1984). We assume the energy-momentum tensor for matter of an ideal gas described as 

T* = p h g u a u p + P g9 ap , (1) 

where u a , po, h g and P g are the four- velocity, the rest-mass density, the relativistic specific enthalpy and the pressure of the 
gas, respectively. The specific enthalpy of the gas, h g , is calculated as h g = (e g + P g )/po where £ g is the energy density of 
the gas. Here, the fluid quantities h g , e g , P g and po are all being measured in the comoving frame of the fluid. The radiation 
stress-energy tensor is calculated from the specific intensity I(x a ; n, v) as 

^ &= 1 J I ^ n ' v "> n °' n0d ' jdn ' ( 2 ) 

where p a = (hv)(l,n) are the photon (or neutrino) four-momentum and n a = p a /hv. Here, the specific intensity I(x a ;n, v) 
is defined for photons (or neutrinos) moving in direction n with the frequency v. The particle number conservation equation 
in the absence of particle creation and annihilation and the conservation equation for the total energy-momentum of gas plus 
radiation are given as 

V a (nu a ) = 0, (3) 
V p {T af) + R afi ) = 0, (4) 

respectively. Here n is the particle number density measured in the comoving frame, and n is related to po as n = po/m g 
where m g is the mass of the gas particle. On the other hand, the radiation four-force density acting on the matter is given as 

G a = ~ / J [*„ J(n, v) - VlJ ]n a dudn, (5) 
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where Xv an d rj u are the opacity and the emissivity, respectively. The invariant emissivity and invariant opacity are r] v /v 2 
and vx v , respectively. The dynamical equations for the matter and the radiation field are described as 

V T al3 = G a , (6) 
Vf)R af) = -G a . (7) 



3 RADIATION HYDRODYNAMICS DESCRIBED BY KERR-SCHILD COORDINATE 

In this study, we assume the background geometry around the rotating black hole written by the Kerr-Schild coordinate where 
there is no coordinate singularity at the event horizon. By using this coordinate, the metric around a rotating black hole is 
described as 

ds 2 = g a pdx a dx 13 , 

= -a 2 dt 2 + yaidx* + f3 % dt)(dx j + f3 J dt), (8) 

where i, j = r, 9, <f> and the nonzero components of the lapse function a, the shift vector /3 l and the metric in three-dimensional 
spatial hypersurface 7^ are given in the geometric unit as 

/ 2mr\' 1 / 2 „_ 2mr/E 
01 = ) '^=1 + 2^7=' 

2mr ^ Asm 2 6 . 2 n (\ 2mr\ 

7rr = 1 + —^-,-yee = E, 7w = — f, — >7r<#. = 74>r = -asm 9 I 1 + I . (9) 

Here, we use the geometric mass m = GM/c 2 , S = r 2 + a 2 cos 2 0, A = (r 2 + a 2 ) 2 - a 2 A sin 2 9 = EA + 2mr(r 2 + a 2 ) = 
E + a sin #(E + 2mr), and A = r — 2Mr + a , where M is the black hole mass, G is the gravitational constant and c is 
the speed of light. The position of the outer and inner horizon, r±, is calculated from A = as r± = m ± (m'-a 2 ) 1 ' 2 . The 
angular velocity of the frame dragging due to the black hole's rotation is calculated as u> — —gt<j>/g4>4, = 2mar/A. 



3.1 Reference Frames and Orthonormal Tetrads 

Three reference frames are used in this study: (1) the Kerr-Schild frame (KSF) which is the frame based on the Kerr-Schild 
coordinate describing the metric, (2) the locally non-rotating reference frame (LNRF) which is the orthonormal frame fixed 
to the coordinate, and (3) the comoving frame where the element of the fluid is at rest. The interactions between the matter 
and the radiation are introduced in the comoving frame. The LNRF is calculated from a stationary congruence formed by 
observers wit h a future-directed unit vector orthogonal to t =constant. The components of the four velocity for such observer 
are given as ( Frolov fc Novikov 199fih 

ti M = — q5* , and it* = a -1 , u = — q -1 /? 1 , (i = r, 9, </>), (10) 

respectively. Since the vorticity tensor vanishes for this congruences, the reference frame formed by this congruence is locally 
non-rotating. Then, the tetrad vectors, eJJ and for the LNRF are given as 

(11) 
(12) 

(13) 

(7^) 1/2 ] , (14) 

(15) 
(16) 
(17) 
(18) 

Here, the hat denotes the physical quantities measured in the LNRF. The base ea = d/dx a of the LNRF can be expressed 
by the coordinate base d/dx 11 as 

L = l(L-gr±\ ± = ^(±.Jr*£\ A = A = ( 19) 

dt a \m P drj> df vr \dr 7w dcj> ) ' d § 09' ^f^dcp' K ' 



t 


[a, 0, 0, 0] , 


r 

e M = 


[wr 1/2 , { YT 1/2 , 0, 0] 


§ 


[0, 0, (yee) 1/2 , 0] , 


4> 

4 = 


[/3 r 7r^»(7^)~ 1/2 , 7r<#,(7^)~ 1/2 


and 




< = 


[o- 1 , - fa' 1 , 0, 0] , 


4 = 


[0, ( 7 rr ) 1/2 , 0, 7 P *(7 PP )" 1/3 ] 




[0, 0, (/ 9 ) 1/2 , 0] , 




[0, 0, 0, ( 7 ^)~ 1/2 ] • 
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where fi = t, r, 9, <j>. Here, (x°, x 1 , x 2 , x 3 ) = (i, f, 9, 0). The components of the fluid's three velocity measured by a 
fiducial observer who is fixed with respect to the coordinates in the LNRF are calculated as 



(i = r, 



(20) 



where u a is the four-velocity in the LNRF. The components of the three velocity are explicitly calculated as 



u 



n + 



7r<£ 



(21) 



where = u^/u* is the angular velocity and we have used u* = — u- t — au*. The Lorentz factor 7 for this three velocity is 
calculated as 



7 = (1 - v ) 



~2n-1/2 t 

= au , 



(22) 



where v = v • v = «;-ti l = t) r + t) e + 0^ . 

A tetrad base for the comoving frame d/dx a is calculated by the Lorentz transformation as 

^ =A * a (v)-^, (23) 
dx a dx? 

where the bar denotes the physical quantities measured in the comoving frame. The components of the Lorentz transformation 
A^(v) are given as A* f = 7, A* f = jv\ A f j = 7U, and A'j = 8 l j + v'vjj 2 /(l + 7) (i, j = r, 9, <t>). Here, (a; 5 , x 1 , x 2 , a: 3 ) = 
(t, f, 6, 0). The components of the base of the comoving tetrad d/dx a can be expressed by the coordinate base d/dx a as 



7 9 



dt 

d_ 

df 



a dt 

7"y r 9 
a dt 



+ 



/T\ _9_ JVg_d_ _ [ 

a J dr . /ZZ 89 



7r0 j _9_ 

7^0 / 90' 



7 rr 1 + 



VrV^ 2 



7 + 1 



-2-2 

«r7 

7 + 1 



Q 



9 


7«e 


9 


90 


Q 


dt 


9 


7«0 


9 


d4> 


a 


dt 



+ 7*e 
+ 1V4, 



J<t>4> 
£V7 \ P 



1 + 



— + — 

dr 



1 + 7 



7ee 

_9_ 

90' 



VrVg'y 2 \ d 



7 + 1 

«r7 \ 



_9 _L 

dr 



1 + 



7+1 y 96» 



*e7 2 \ d vgj 



7 + 1 



+ - 



7 + 1 



■ «rV7 



7r0 \ d_ 

7<W ) d< t>' 



7 + 1 I a 



d 1 / ue«07 2 ^ 9 
9 r ^7+1 



1 + 



~2 ~ 2 

^7 



7r0 / VrVtpJ 2 



14>4> \ 7 + 1 



d 



(24) 



In the similar manner, the inverse transformation from the base of the comoving tetrad to the coordinate base is calculated by 
using the inverse Lorentz transformation A a ^(— v), and the coordinate base d/dx a is calculated from the base of the comoving 
tetrad d/dx a as 



ld_ 

a dt 



= 7 



dr 



1 d 
P^~ 6 d9 



1 - 



«07r0 



70<* 



d_ 

dt 



-^Vr + 



fir 



7 2 t) 2 

1+-? 



Vr 7 V 7 + iy qVt^ V 7 



'+1 



_9 

9f 



-7«e + 



-7#* + 



/3 r / 7 2 u r u£)\ /3 r 7r^ / 7 "e«0 



Vr 7 V 7 + 1 / V 7 + 



f3 r ( y 2 v r v<f, \ f3 r -y r <p 



1 + 



-2*2 

7 ^ 



a^fr\l + 1 J ay^\ 7+1 

-7 I Or + «07i 



_9_ 



7 rr \ d 



14>4> J dt 



■ + 



1 + 



-2-2 

7 «r 

7 + 1 



+ 7r^>1 



7" / y 2 v r v$ 



j 2 v r ve 
7 + 1 



+ 7r 




T+Ty 1 



^ + 



y 2 VrV<f> 

7 + 1 



\ 7+1 

+ 7 



_9 




1 + 



7 2 ^ 
7 + 1 



90' 



~- d f 7 2 «- 



0r0g\ _9_ / 7 2 «a \ _^ / 7 2 «t>^ \ _9^ 
+ 1 J df \ 7 + V 9(9 V 7 + 1 y 90' 
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l d 



8 ( ^VrV,!, 



700 



'di 



7 + 1 



d_ 

Or 



7 + 1 



^ + 1 



-2-2 

7 



d 



(25) 



3.2 Radiation moments 

The radiation energy density E, the radiation flux F l and the radiation pressure tensor P %J are defined as the zeroth, the first 
and the second moments of the specific intensity I v (x^ , n), respectively. We denote the radiation moments as 



E = / / IodOdn, F 1 = / / hndvdQ., P lJ = 
when measured in the LNRF, and 
E - I I hdvdQ., F l = I I hffdvdn, P ij = 



Infin 3 dvdQ,, 



Iunn 3 dvdQ,, 



(26) 



(27) 



when measured in the comoving frame. Correspondingly, the radiation stress tensors for the LNRF and the comoving frame 
are given as 



R 



( k 


F r 


pe 


F^ \ 


F r 


prr 


pre 


pr<t> 


F 6 


pr6 


pee 


p»0 


\ P* 


pr<t> 


pe<t> 


p4><t> j 



and R 



a/3 



f 


E 


F r 


po 


F* \ 




F r 


prr 


pre 


pr<j> 




F e 


pre 


pee 


pe<t> 




p<!> 


pr<t> 


pe<t> 


p4>4> j 



(28) 



respectively. The contravariant components of the radiation stress tensor R a/S are calculated from R al3 by the transformation 
as 

p a P ^ x pp,0 

dx^ dx" 

and explicitly given as 



(29) 



R° 



E 



-yrr p r 



-E 



Y r F r 



1 

Q 



-E 



OL\/ ^jee 
v r . p4> 



^rr prr 



2f3 r 



,yrr pr ^_ 



E 



F r + 





1 (- 




\fyee ^ 


pr<t> 


h 7 A 


\Jl4>4> 


\/ r y rr 



ryrr pre _ M_ pt 

a 



prr _ r_ 



F r + 



p4 



F 



ay -fee 



a \ 



P \ 



Sip* 



jrd> 



7 



rtf> 



pee 
7ee 



+ 

7^0 V7' 

1 / P e0 



700 / 

7' 



a 



+ -+=P r 



t4> 



F r + 



p<t> \ 



7<W J 



pS4> 

+ -+=P r 



+ 



2rf* hr<j> , ( 7 ^) 2 



700 y/ r y rr ~. 



__pr<P + 



r 



(30) 



By using the Lorentz transformations, the radiation moments measured in the comoving frame are calculated from those 
measured in the LNRF as 



R 



a/3 



and explicitly given as (jMihalas fc Mihalaslll98l |Parkll2006h 



F i = -7V.E + 7 



SU 7 + 



7 + 1 



(31) 

(32) 
(33) 
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Yvv J E - 7 v l 5l + v>dl + 2^—v % v ] v k ) F k + [51 + - 



v'Vk I I S 3 , 

7 + 1 I V 7+1 / V 7+ 1 



V J t>! P 



(34) 



Inversely, the radiation moments measured in the LNRF are calculated from those measured in the comoving frame as 
R & ? =A%(v)Ai(v)R^, (35) 
and explicitly given as 



E = 7 (E + 2v i F t + v l v J P lJ ) , 



F* = j 2 v 4 E + j 



5*+ 7 + 



7 + 1 



7+1 / V 7+1 



5! 



7 + 1 



P 



(36) 
(37) 

(38) 



3.3 Radiation four-force density 



The radiation four-force density measured in the comoving frame is given as ( Mihalas fc Mih alas 1984) 



G a = - 

c 



(39) 



where \ an d V are the mean opacity and the emissivity measured in the comoving frame, respectively. The time component 
G has the dimension c _1 times the net rate of the radiation energy per unit volume, and the spatial component G % has the 
dimension of the net rate of the momentum exchange between the matter and t he radiation. The time component o f the 
radiation four- force density measured in the comoving frame can be calculated as (Mi halas fc Mihalas 1984 ; Park 20061) 



G = T - A 

where the heating function f and the cooling function A are defined as 



1 



rjdvdCl. 



The components of the radiation force G°" is calculated from those in the comoving frame G a by the transformation 

Q a G^ 

and explicitly given as 

G* = l((? + ViG 
a \ 

\ff^G f + 7 
1 



(40) 



(41) 



(42) 



G r 
G ti 



7 V r 



^)G r +7 



7" ——- r \v l G , 

7 + 1 a 



G" = 



\Jloe 



r<e , - ~ rH i 7 2 «S - rA 

G + 7« e G + -ViG 

7 + 1 



+ 



G p + 7 



^ , 7 



G* + 



7 + 1 



fi.G 1 



(43) 



4 RADIATION HYDRODYNAMIC EQUATIONS 
4.1 Continuity equation 

As a continuity equation, now we consider the particle number conservation, (nu a )- a — 0. This equation can be calculated as 
(ay^nu^j + (ay^m/) = 0, (44) 
where 7 = det7;,j. In the case of the Kerr metric written by the Boyer-Lindquist coordinate 7 = E(E + 2mr) sin 2 8. 
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4.2 Hydrodynamic equations 

The relativistic Euler equations are obtained by the projection of the equation of the energy momentum conservation V pT al3 — 
G a on the specific directions by using the projection tensor P a/3 = g a/3 + u a vP as Pp\7 ' sT 0S — PpG 13 . From this, we can 
obtain p o h g u V u a + {g af) + u a u )d,3P g = G a + u a upG 9 . 
The Euler equations in r, 9 and (^-directions are given as 

pohgU^- + PohgU^- + pahgY^u^u' + (g tr + u r u)^- + (g" + u T v?)-jpr = G r + v^upG 13 , (45) 

MX^ + MX^ + Pofc^u^+ttV^ + ^+iiV)^ = G e +u e u G , (46) 

p ohg u t ^+p h g u^+p h g r^yu'' + (g^ + u r u t )^ + (g^+u r u^ = G* + W %G", (47) 
where 

Fp^u' = r" t (u) 2 + r" r (u) 2 + Yg e (u) 2 + r^,(?/) 2 + 

+2r t ™ wV + 2rg,uV + 2I%uV + 2T? e u r u e + 2r^uV + 2r^«V> (48) 

and now a = r, 8, <j>. The Christoffel symbols are given in App.fX] u a G a is calculated by the heating and cooling function 
defined in the comoving frame as 

u a G a = ~G l = A - r. (49) 
This is also calculated as 

u a G a = [(-a + jS^JG* + jr4,p T [G T + G 4 ")] u* + [irr{G r + /TG 4 ) + lr4> G 4 '] u r + leB G u e 

+ [f^G* + lr4 ,(G r + /3 r G 4 )] u*. (50) 
The he local energy conservation is obtained from u a X7 pT a0 — u a G a calculated as 



4.3 Radiation moment equations 

The radiation moment equation \7pR alB = —G a gives the equation for the energy density, the radiation flux and the radiation 
pressure tensor. The radiation energy equation is obtained from f-component of this equation V ' a R ta = — G* calculated as 



dR u + 1 
dt ct\Jl 



d_ 

dr 



+ 



-G l 



The radiation momentum equation in r-direction X7 a R ra = — G r is calculated as 



dR tr 1 
+ 



d_ 

dr 



+ 



OR 



pi" TtP-f _ sir 



dt ' cty^y 

In the similar manner, the radiation momentum equation in 0-direction V ' a R 0a = — G e is calculated as 



dR w 1 
+ 



dt 



|(«^**)+|(«v-«- 



dR H 
d<f> 



r e R /5 7 _ r ,o 

1 PySt lJT 



The radiation momentum equation in ^-direction V a R = — G^ is calculated as 



dR^ 1 
+ 



dt 



d 



+ 



dcj> 



4> T)Pl 



-G* 



(52) 



(53) 



(54) 



(55) 



Finally, by inserting the components of the radiation stress tensor given by Eq. l|30j> into Eqs. J52J, (|53|l , (|54|l and ((55} , we 
obtain the radiation moment equations as 



d_( E 

dt 



+r 



i d 

oty/^y dr 



+ 2FI- 

a 



a 



Q^F r - —E) +2Y\ e 
a 



+ 



1 d 

t^f de 



d 

+ TTT 



7 



r<j> 



F r 



P0 



F 



7 r T r - ^- yfrF r + ( — ) E 



a y ■fee 
poo 

J 00 



p4>4> 



., r 4 > 



F r + 



p<t> \ 



+ 



27'' .p^ + (7 r0 ) 2 



7^0 \/ Y T l<t>4> 
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"fee ^ Q 



+ 2r 



rep 



pr0 ^r<p 



+ 

700 \ll 



P ( y0 ^ + ^ 



'700 



+2r£ 



6cp 



p 9 <t> Y<t> „ 



9i 



Tee \ V 700 



= -G 4 , 



(56) 



a V a / 



+ 



1 9 



+ 



1 9 
+1* 



a 



«V7 



a 

pr<j> 



1H> V7 



F 

/3 r / 7 r * 



F r + 



+ 2r[ r 



yrprr _ W_ /777- ' pr + / ^\ 

a V a / 



+ 2rl 



F 



+ 2r 



\ ay 7ee 
pee 
"fee 



■ 00 



tip 



1 7 



F r + 



F^ \ 



+ 



+2r. 



lee 
1 



+ 2r 



60 



_9 [ _F 

9t 



7ee \ V 



ay "lee 



rep 



= -G r 
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where the ChristofFel symbols are given in App. [X] 



5 CONCLUDING REMARKS 

In this study, we have derived equations of fully general relativistic radiation hydrodynamics around a rotating black hole 
by using the Kerr-Schild coordinate where there is no coordinate singularity at the event horizon. Both the matter and the 
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radiation are affected by the frame-dragging effects due to the black hole's rotation. Since the radiation usually interact with 
matter moving at relativistic velocities near the event horizon, the interaction between the radiation and the matter should be 
treated fully relativistically. This can be done if the interplay between the radiation and the matter is evaluated in the comoving 
frame where the matter is at rest. In this approach while the interactions between matter and radiation are introduced in the 
comoving frame, the equations and the derivatives for the description of the global evolution of both matter and the radiation 
are given in the Kerr-Schild frame (KSF) which is a frame fixed to the coordinate describing the central black hole. As a frame 
fixed to the coordinate, we use the locally non-rotating reference frame (LNRF) which is a radially falling frame when the 
Kerr-Schild coordinate is used. We can see in the derived equations that both the matter and the radiation are affected by the 
frame-dragging effects due to the black hole's rotation. It is widely known that the moment equations truncated at the finite 
order do not produce the complete system of equations, i.e., the number of v ariables are larger than the number of equations. 
So, the additional equations are required to close the system of equa tions ( Chandrasekharl 1960l ; Mihalasl 197p| ; Pomraning 
19731 ; iRvbicki fc Lightmanlll979l ; llVIihalas fc Mihalasl 1984 ; IShul 119911 ). These additional equations should be also described 



fully relativistically when the radiation interact with the matter moving at relativistic velo cities. As one possible way with 
respect to such additional equations, the flux-limited diffusion approximation proposed by iLevermore fc Pomraning] (|1981f ) 
is sometimes used. The covariant theories for the flux-limited diffusion ap proximation is presented for radiation propagating 
through inhomogeneous and nonstationary media I Anile fc Romanol 1992h . It is noted that constructing the closure relations 
which are essentially required for the radiation hydrodynamical calculations is non-trivial when the photons and the matters 
are relatively relativistically moving. This is mainly because the truncation of the moment equations at the finite order can 
not be performed frame-independently. So, in such cases, the equations given in this articles are of no use if the sufficient 
closure relations can not be constructed. In the case of the calculatio ns based on the K err-Schild coordinate, the fluid velocity 
is usually smaller than the speed of light even at the event horizon ([Takahashi 2007a), for example, the gamma factor of the 
fluid at the horizon 7 < 1-100. So, if the photon motion have the speed similar to the fluid's speed and the closure relations 
which is correct for such motions can be constructed, the equations given in this study can be used for the description of such 
radiation hydrodynamic flows. It is also note that when effects of the radiation stress can not be neglected, the dissipation 
effects can not be correctly treated so far. This is mainly because we do not know the covariant theory describing the local 
dissipation effects such as viscosity. One of the pr omising approach is g iven by the method based on the extended causal 
thermodynamics such as the Israel-Stewart theory ( Israel fc Stewartlll979l ) where the causality violating infinite signal speeds 
are eliminated (see also, Anile, Pavon fc Romano 1998). However, while the hydrodynamical equations based on such theory 
are formulated by iPeitz fc App] l|l998l) . there is no formulations based on such theory for the radiation hydrodynamic or 
magnetohydrodynamic equations. This is one of the important studies in terms of the relativistic radiation hydrodynamics in 
future. 
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APPENDIX A: DIFFERENTIAL VALUES OF METRIC COMPONENTS AND CHRISTOFFEL 
SYMBOLS 

Here, we present the explicit expressions for the metric components and the their differential values with respect to r and 9 
used in this paper. Non-zero components of the metric are given as 
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The differential values of the metric components are given as 
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where we have used 

9 r E = 2r, d r A = 2 [2rE + (r + m)a 2 sin 2 6»] , 9 e E = -a 2 sin 29, d g A = -a 2 A sin 29. 
The Christoffel symbols (a = t, r, 9, <f>) are given as 
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